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ABSTRACT 

We  present  spectroscopy,  astrometry,  and  photometry  of  the  brown  dwarf  WISEP  J004701. 06+680352.1 
(W0047+68),  an  unusually  red  field  L  dwarf  at  a  distance  of  12.2  =b  0.4  pc.  The  three-dimensional  space  mo¬ 
tion  identifies  it  as  a  member  of  the  AB  Dor  Moving  Group,  an  identification  supported  by  our  classification  of 
W0047+68  as  intermediate  surface  gravity  (INT-G)  using  the  Allers  &  Liu  near-infrared  classification  system.  This 
moving  group  membership  implies  near-solar  metallicity,  age  ~  100-125  Myr,  M  ~  0.018  M0,  and  logg  ~  4.5; 
the  thick  condensate  clouds  needed  to  explain  the  infrared  spectrum  are,  therefore,  a  result  of  surface  gravity  that 
is  lower  than  that  of  ordinary  field  brown  dwarfs.  From  the  observed  luminosity  and  evolutionary  model  radius,  we 
find  reff  ~  1300  K,  a  temperature  normally  associated  with  early  T  dwarfs.  Thick  clouds  are  also  used  to  explain 
the  spectral  properties  of  directly  imaged  giant  planets,  and  we  discuss  the  successes  and  challenges  for  such 
substellar  models  in  matching  the  observed  optical  and  infrared  spectra.  W0047+68  shows  that  cloud  thickness  is 
more  sensitive  to  intermediate  surface  gravity  than  in  most  models.  We  also  present  a  trigonometric  parallax  of  the 
dusty  L6  dwarf  2MASS  J2148 1628+4003593.  It  lies  at  8.060  zb  0.036  parsecs;  its  astrometry  is  consistent  with  the 
view  that  it  is  older  and  metal-rich. 

Key  words:  brown  dwarfs  -  infrared:  stars  -  stars:  individual  (WISEP  J004701. 06+680352.1, 

2MASS  J21481628+4003593) 


1.  INTRODUCTION 

One  of  the  key  challenges  in  substellar  astronomy  is  to  use 
observations  of  colors,  absolute  magnitudes,  and  spectra  to 
deduce  the  fundamental  physical  parameters  of  brown  dwarfs 
and  gas  giant  planets.  Determining  the  luminosity,  effective 
temperature,  surface  gravity,  composition,  mass,  and  radius 
would  allow  many  astrophysical  applications,  including  testing 
formation  scenarios  for  “planetary  mass  objects,”  but  is  difficult 
in  part  due  to  the  importance  of  condensate  cloud  formation 
at  substellar  atmospheric  temperatures.  The  vast  majority  of 
the  nearly  two  thousand  known  field  brown  dwarfs  can  be 
described  by  a  well-defined,  one-dimensional  spectral  sequence 
of  late-M,  -L,  -T,  and  -Y  dwarfs,  with  a  dramatic  change  in 
cloud  properties  at  the  L/T  transition  around  1400  K  (see  the 
review  of  Kirkpatrick  2005).  It  is  now  clear,  however,  that 
many  of  the  directly  imaged  planetary  mass  companions  have 
significant  differences  from  ordinary  field  brown  dwarfs.  The 
companions  2M  1207b  (Chauvin  et  al.  2004;  Patience  et  al. 
2010)  and  HR  8799b,  c,  d,  and  e  (Marois  et  al.  2008,  2010), 
are  considerably  redder  than  field  brown  dwarfs  of  similar 
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temperature  or  luminosity.  The  planet  Pic  b  has  apparently 
normal  near-infrared  colors  for  a  mid-L  dwarf  but  very  red  mid- 
infrared  colors  (Currie  et  al.  2013),  and  the  L'  magnitudes  of 
the  HR  8799  planets  are  also  brighter  than  expected  (Skemer 
et  al.  2012).  Spectra  of  the  HR  8799  planets  show  considerable 
diversity,  which  may  require  non-equilibrium  mixing,  different 
cloud  properties,  and  perhaps  enhanced  metals  or  non- solar 
abundances  (Bowler  et  al.  2010;  Barman  et  al.  2011;  Marley 
et  al.  2012;  Konopacky  et  al.  2013;  Oppenheimer  et  al.  2013), 
but  distinguishing  between  these  possibilities  is  difficult. 

Rarer  classes  of  field  brown  dwarfs,  especially  the  unusually 
dusty  ones,  offer  a  way  forward,  as  they  can  be  more  accessible 
to  observations  than  planets.  A  number  of  systematic  studies 
of  young  (<150  Myr)  late-M  and  early-L-type  brown  dwarfs 
(Kirkpatrick  et  al.  2008;  Cruz  et  al.  2009;  Patience  et  al.  2012; 
Faherty  et  al.  2013a,  2013b;  Allers  &  Liu  2013a;  Manjavacas 
et  al.  2014)  have  demonstrated  their  red  colors  and  other 
spectral  peculiarities,  which  can  be  ascribed  to  low-surface- 
gravity  dusty  atmospheres.  Unfortunately,  even  cooler  brown 
dwarfs  are  necessarily  fainter  and  more  difficult  to  discover. 
Gizis  et  al.  (2012,  hereafter  Paper  I)  noted  that  there  were  then 
only  five  known  field  L  dwarfs  with  J  —  K  >  2.3,  and  none  had 
trigonometric  parallaxes.  The  field  is  developing  rapidly  with 
new  sky  surveys,  and  the  two  currently  reddest  (/  —  K  =  2.8) 
field  L  dwarfs  known  are  PSO  J318.5338-22.8603  (Liu  et  al. 
2013b)  and  ULAS  J2227 11-004547  (Marocco  et  al.  2014).11 


11  All  names  are  abbreviated  after  first  use:  PSO  J3 18-22  for  PSO 
J318.5338-22.8603,  W0047+68  for  WISEP  J004701. 06+680352,  etc. 
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Figure  1.  Moderate-resolution  near-infrared  IRTF  spectrum  of  W0047+68  compared  to  other  unusually  red  L  dwarfs  from  Allers  &  Liu  (2013a). 


PSO  J318-22’s  trigonometric  parallax,  tangential  velocity  and 
unusual  spectrum  suggests  it  is  a  planetary-mass  member  of 
the  ft  Pic  Moving  Group  (Zuckerman  et  al.  2001),  and  in  many 
ways  is  similar  to  2M  1207b.  It  has  the  faintest  J- band  absolute 
magnitude  (14.8)  of  any  L  dwarf  other  than  2M1207b  (16.13, 
Gizis  et  al.  2007).  Besides  PSO  J3 18-22,  the  best  studied  of 
the  extremely  red  L  dwarfs  are  2MASS  J2 148 1628+4003593 
(Looper  et  al.  2008),  which  is  thought  to  be  old  (>200  Myr), 
2MASSW  J22443 16+204343  (Dahn  et  al.  2002),  which  has 
been  proposed  to  be  young  (Kirkpatrick  et  al.  2008),  and 
2MASS  J035523. 37+1 13343.7,  a  dusty,  red  L-type  dwarf  with 
low  surface  gravity  (Cruz  et  al.  2009;  Faherty  et  al.  2013b;  Liu 
et  al.  2013a). 

WISEP  J004701. 06+680352  (Paper  I)  was  the  reddest  (/  - 
Ks  =  2.55  it  0.08)  known  L  dwarf  at  the  time  of  discovery. 
Although  its  low-resolution  near-infrared  spectrum  resembles 
older  L  dwarfs  (see  Paper  I),  Thompson  et  al.  (2013)  have 
suggested  that  it  is  young  and  has  low  surface  gravity  on  the 
basis  of  its  moderate-resolution  spectrum,  and  Gagne  et  al. 
(2014)  argue  that  it  and  2M2244+20  are  probable  members  of 
the  AB  Dor  Moving  Group  (ABDMG).  In  this  paper,  we  present 
new  observations  and  analysis  of  the  brown  dwarf  W0047+68, 
including  optical  and  near-infrared  spectra,  photometry,  and  a 
preliminary  trigonometric  parallax,  as  well  as  a  parallax  for 
2M2 148+40.  We  argue  that  W0047+68  is  a  bona  fide  member  of 
the  AB  Dor  Moving  Group  (ABDMG),  and  confront  theoretical 
dusty  models  with  our  new  observational  constraints.  We  find 
that  2M2148+40’s  space  motion  is  consistent  with  the  view  that 
it  is  older  and  metal-rich. 

2.  DATA  AND  OBSERVATIONS 
2.7.  Near-infrared  Spectra 

A  moderate-resolution  near-infrared  spectrum  covering  the 
range  of  0.8  to  2.4  /xm  was  obtained  using  the  IRTF  SpeX 


spectrograph  (Rayner  et  al.  2003)  in  cross-dispersed  mode  (here¬ 
inafter  SXD)  on  UT  Date  2012  September  26.  The  resolution  is 
1200.  The  data  were  reduced  using  the  facility  pipeline  Spextool 
(Cushing  et  al.  2004).  Telluric  features  were  corrected  using  a 
nearby  A  star  using  the  method  of  Vacca  et  al.  (2003).  Spec¬ 
tra  with  the  same  observing  setup  were  used  by  Allers  &  Liu 
(2013a)  to  develop  a  spectroscopic  classification  system  that  is 
sensitive  to  surface  gravity:  the  key  indicators  include  the  Ki 
lines  in  the  J  band  and  the  FeH  bands  in  J  and  H  bands.  In 
Figure  1,  we  compare  our  spectrum  to  other  unusually  red  L 
dwarfs  2M2244+20, 2MASS  J01033203+1935361  (Kirkpatrick 
et  al.  2000;  Faherty  et  al.  2012),  and  2M2148+40.  W0047+68 
is  similar  to  2M2244+20,  but  has  considerably  lower  surface 
gravity  than  2M2148+40.  Applying  the  Allers  &  Liu  (2013a) 
system,  we  classify  W0047+68  as  L7  INT-G  (Table  1). 

High-resolution  near-infrared  spectra  of  W0047+68  were 
obtained  with  Keck/NIRSPEC  (McLean  et  al.  1998,  2000)  on 
UT  Date  2013  September  17  to  determine  W0047+68’s  radial 
velocity  and  projected  rotational  velocity  (vsini).  Data  were 
obtained  using  the  N7  filter  and  (L432  x  12"  slit  (3  pixel 
dispersion),  providing  1.99-2.39  /xm  spectroscopy  in  seven 
orders  at  a  measured  resolution  of  25,000  (Av  =  12kms_1). 
Two  dithered  exposures  of  1200  s  were  obtained  at  an  airmass 
of  1.50,  followed  by  observations  of  a  nearby  A  star.  Data 
in  the  order  of  33,  spanning  2.293-2.319  /xm,  were  extracted 
and  analyzed  due  to  the  strong  CO  bands  present  in  this 
band.  We  used  a  combination  of  REDSPEC12  and  custom 
IDL  routines  that  perform  a  forward-modeling  Monte  Carlo 
Markov  Chain  analysis  of  the  extracted  spectrum  using  solar- 
metallicity  BT-Settl  models  (Allard  et  al.  2011)  and  a  telluric 
absorption  spectrum  from  Livingston  &  Wallace  (1991);  this 
analysis  will  be  described  in  detail  in  a  forthcoming  paper  (A.  J. 
Burgasser  et  al.,  in  preparation;  see  Blake  et  al.  2010  for  a 
similar  algorithm).  Figure  2  displays  our  extracted  spectrum 
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Figure  2.  Analysis  of  our  Keck/NIRSPEC  spectrum  of  W0047+68  (black  line).  Left:  the  best-fit  model  (red  line;  Teff  =  2300 K  and  logg  =  5.5  cgs)  has  been 
broadened  by  v  sin  i  =  2  km  s-1  and  shifted  by  a  heliocentric  radial  velocity  of  —21 .5  km  s_1 .  The  green  line  shows  the  combination  of  model  and  telluric  absorption 
(from  Livingston  &  Wallace  1991)  which  best  fits  the  data.  Residuals  are  shown  at  bottom  (black  line)  and  are  comparable  to  the  noise  (±l<x  indicated  by  the  gray 
region  around  zero  flux).  Right:  a  comparison  of  two  BT-Settl  (Allard  et  al.  2011)  models,  Teff  =  2300  K  with  logg  =  5.5  and  Teff  =  1400 K  with  logg  =  4.0, 
illustrating  that  high-resolution  spectra  in  this  wavelength  region  can  be  used  for  velocity  information  but  do  not  reliably  constrain  temperature. 


Table  1 

WISEP  J004701. 06+680352.1 


Parameter 

W0047+6803 

Remarks 

Sp.  type 

L7  INT-G 

Infrared,  Allers  &  Liu  (2013a)  system 

Sp.  type 

L7  pec 

Optical,  Kirkpatrick  et  al.  (1999)  system 

^abs 

82  ±3 

mas 

prel 

434  ±3 

masyr-1 

e 

117.0  ±0.5 

degrees,  east  of  north 

brad 

-20.0  ±  1.4 

kms-1 

v  sin  i 

4.3  ±  2.2 

kms-1 

z 

18.95  ±0.03 

AB  mag,  Pan-STARRS 

y 

18.12  ±0.03 

AB  mag,  Pan-STARRS 

u 

-8.6  ±  1.0 

kms-1 

V 

-27.8  ±  1.3 

kms-1 

w 

-13.6  ±0.5 

kms-1 

Mj 

15.17  ±  0.11 

2MASS 

Mh 

13.54  ±  0.09 

2MASS 

MKs 

12.62  ±  0.09 

2MASS 

MWi 

11.47  ±0.08 

WISE 

Mw  2 

10.82  ±  0.08 

WISE 

MW3 

9.71  ±0.10 

WISE 

mbo\ 

16.30  ±  0.03 

Section  2.5 

L/Lq 

(3.58  ±  0.29)  x  10-5 

Section  2.5 

BCj 

0.70  ±  0.08 

BCk 

3.25  ±  0.04 

Mj  (MKO) 

15.05  ±0.11 

Color  term  from  spectrum,  see  Paper  I 

Mh  (MKO) 

13.65  ±  0.09 

See  Paper  I 

Mk  (MKO) 

12.56  ±  0.09 

See  Paper  I 

MM' 

10.6 

Using  2M2244,  see  Section  2.4 

If  AB  Dor  Association  member 

Teff 

-1270 

K,  see  Section  3.3 

Age 

-120 

Myr 

Mass 

-0.018 

Mq 

log# 

-4.5 

If  older  Galactic  Disk  member 

Teff 

-1410 

K;  see  Section  3.3 

Age 

^1000 

Myr 

Mass 

<0.055 

Mq 

logg 

-5.0 
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Wavelength  (microns) 

Figure  3.  Optical  spectra  of  W0047+68,  along  with  the  red  L6.5  (opt)  2M2244+20  from  Kirkpatrick  et  al.  (2008)  and  the  red  L6  (opt)  2M2 148+40  from  Looper  et  al. 
(2008).  All  four  spectra  have  been  smoothed. 


and  the  best-fit  model.  The  individual  best-fit  model  has  nrad  = 
—21.5  km  s_1and  v  sin  i  =  2  km  s_1,  but  marginalizing  over  the 
fit  parameters,  including  other  temperatures  and  gravities,  we 
find  v  sin  i  =  4.3  zb  2.2  km  s-1and  a  heliocentric  radial  velocity 
of  —20.0  ±1.4kms_1,  which  we  adopt  as  our  best  estimates 
for  W0047+68.  Our  modeling  requires  higher  resolution  than  is 
publicly  available  for  the  full  BT-Settl  grid;  the  subset  of  the  grid 
that  we  use  has  Tes  >  2000  K,  hotter  than  we  believe  W0047+68 
can  possibly  be.  However,  these  CO  lines  are  not  very  sensitive 
to  temperature.  We  show  in  Figure  2(b)  a  comparison  of  the 
Tiff  =  1400  K,  log  g  =  4.0  model  (which  we  discuss  later)  to 
the  hotter  models  used  by  the  modeling  routine  at  the  highest 
resolution  available.  The  lines  are  very  similar,  suggesting  that 
the  radial  velocity  and  v  sin  i ,  but  not  the  temperature  or  gravity, 
can  be  extracted  from  the  Keck  data. 

2.2.  Optical  Spectra 

Optical  (far-red)  long-slit  spectra  were  obtained  with  the 
MMT  and  Gemini-North  telescopes.  The  MMT  observations 
were  on  UT  Date  2012  August  26  with  the  Red  Channel 
spectrograph  with  grating  270,  using  three  600  s  exposures. 
Conditions  were  non-photometric.  This  time  was  allocated  by 
NOAO  as  program  2012B-0233.  The  wavelength  coverage  was 
6170-9810  A  with  a  resolution  of  ~12  A,  but  we  make  no  use 


of  the  spectrum  redder  than  9250  A  due  to  strong  telluric  water 
absorption.  The  Gemini-North  observations  (Gemini  program 
GN-2012B-Q-105)  were  on  UT  Date  2012  September  11  with 
the  GMOS  spectrograph  (Hook  et  al.  2004)  using  grating  R831. 
Four  600  s  exposures  were  taken  during  cloudy  conditions.  The 
wavelength  coverage  was  6340-8460  A  with  a  resolution  of 
~2A,  but  for  all  plots,  we  smooth  to  a  resolution  of  ~10A 
due  to  low  signal-to-noise.  The  spectra  were  processed  using 
standard  IRAF  tasks. 

The  MMT  and  Gemini  spectra  are  shown  in  Figure  3,  along 
with  2M2244+20,  classified  as  L6.5  in  the  optical  (Kirkpatrick 
et  al.  2008),  and  2M2148+40,  classified  as  L7  in  the  optical 
(Looper  et  al.  2008).  The  core  of  the  Ki  doublet  and  the  CrH 
and  FeH  features  are  shown  in  Figure  4  along  with  the  original 
Kirkpatrick  et  al.  (1999)  L6-L8  standard  star  spectra.  The  CrH 
and  FeH  features  most  strongly  support  an  L7  classification, 
and  the  Ki  doublet  and  overall  appearance  of  the  spectrum  is 
consistent  with  this.  It  is  evident,  however,  that  the  Rbi  and 
Cs  i  lines  are  much  weaker  in  W0047+68  than  in  the  standard, 
or  typical,  L6-L8  dwarfs  (Kirkpatrick  et  al.  1999,  2000), 
suggesting  that  a  classification  of  “L7  peculiar”  is  appropriate. 
In  Figure  5,  we  show  the  same  comparison  to  the  unusually  red 
L  dwarfs  2M0355+11  (Cruz  et  al.  2009),  2M2148+40  (Looper 
et  al.  2008),  and  2M2244+20  (Kirkpatrick  et  al.  2008).  These 
spectra  were  classified  as  “L5y,”  L6,  and  L6.5  respectively. 
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Figure  4.  Optical  MMT  spectrum  of  W0047+68  (black)  compared  to  the  Kirkpatrick  et  al.  (1999)  spectra  of  L6,  L7,  and  L8  standards  in  the  Ki  doublet  region  and  the 
CrH/FeH  region;  on  the  left,  the  spectra  are  normalized  by  the  mean  value  in  the  range  of  0.72-0.81  /zm  and  on  the  right,  0.855-0.860  /zm.  W0047+68  best  resembles 
the  L7  standard,  but  the  Rb  i  and  Cs  i  lines  are  both  weaker  than  expected.  Note  that  the  W0047+68  MMT  spectrum  shown  here  was  taken  with  a  different  telescope 
and  instrument,  but  we  have  smoothed  the  data  to  match  resolutions. 


The  weak  Csi  is  also  evident  in  2M2148+40  and  2M2244+20. 
W0047+68’s  alkali  lines  are  stronger  than  in  the  low  gravity 
2M0355+1 1,  but  2M0355+1 1  is  a  warmer  type.  Given  our  near- 
infrared  classification  of  INT-G,  the  atomic  line  peculiarities  hint 
that  W0047+68  might  be  classified  as  L7 in  the  optical,  but 
we  believe  that  extending  the  Cruz  et  al.  (2009)  system  to  such 
a  late-type  object  would  be  premature  without  more  objects. 
Both  W0047+68  observations  show  lithium  absorption,  but  no 
detectable  Ha  emission. 


2.3.  Astrometry 

W0047+68  is  being  observed  as  part  of  the  USNO  optical 
parallax  program,  which  is  described  in  Dahn  et  al.  (2002). 
Based  on  23  frames  taken  in  two  seasons,  the  preliminary 
absolute  parallax  is  82  =b  3  mas  with  relative  proper  motion 
ti  =  434  zb  3masyr_1  at  position  angle  0  =  117.0  ±  0?5. 
We  derive  the  U,  V ,  W  space  velocity  of  W0047+68  using 
the  observed  radial  velocity  and  astrometry  and  the  IDL  code 
gal_uvw.pro ,  except  that  we  adopt  the  traditional  astronomer 
convention  that  U  is  positive  toward  the  Galactic  center.  Note 
that  W  is  strongly  constrained  (±0.5  kms-1)  by  the  astrometry 
because  of  the  position  on  the  sky  (/  =  122.5,  b  =  +05.2).  We 
note  that  the  absolute  magnitude  at  the  H  band  agrees  with  that 
of  typical  L7  dwarfs  (Dupuy  &  Liu  2012);  W0047+68  is  fainter 
(“sub-luminous”)  at  shorter  wavelengths  (/  band)  and  brighter 
(“over-luminous”)  at  longer  wavelengths  ( K ,  WISE). 


Table  2 

2MASS  J2148 1628+4003593 


Parameter 

2M2 148+40 

Remarks 

Sp.  Type 

L6  FLD-G 

Infrared,  Allers  &  Liu  (2013a) 

Sp.  Type 

L6 

Optical,  Looper  et  al.  (2008), 

^abs 

124.07  ±  0.55 

mas 

/Ael 

901.9  +  0.3 

mas 

u 

-34.0  ±  0.5 

kms-1 

mbo\ 

15.12  +  0.03 

Section  3.3 

L/Lq 

(4.64  +  0.14)  x  10-5 

Section  2.5 

Mj 

14.62  ±  0.03 

2MASS 

Mh 

13.25  ±  0.03 

2MASS 

MKs 

12.23  ±  0.03 

2MASS 

Mw  l 

11.21  ±0.02 

WISE 

Mwi 

10.70  ±  0.02 

WISE 

MW3 

10.13  ±0.06 

WISE 

Few  trigonometric  parallaxes  of  unusually  red  L5  or  later 
dwarfs  exist,  complicating  comparisons  of  W0047+68  to  other 
objects.  We  are  able  to  also  report  a  preliminary  USNO  optical 
parallax  of  2M2 148+40,  which  was  classified  as  L6  in  the  optical 
and  L6.5  peculiar  in  the  near-infrared  by  its  discovery  paper 
(Looper  et  al.  2008)  and  is  L6  FLD-G  in  the  Allers  &  Liu  (20 13a) 
system.  The  observed  absolute  parallax  is  124.07  ±0.55  mas, 
with  a  proper  motion  of  901.9  ±  0.3  mas  yr-1  at  an  angle  of 
59?  1  ±  0?  1  (Table  2).  The  distance  of  8.060  ±  0.036  pc  places  it 
just  outside  the  “8  pc”  sample.  The  radial  velocity  is  unknown, 
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Figure  5.  Same  as  Figure  4,  but  for  the  unusually  red  L  dwarfs  2M0355+11  (Cruz  et  al.  2009),  2M2 148+40  (Looper  et  al.  2008)  and  2M2244+20  (Kirkpatrick 
et  al.  2008).  W0047+68  appears  later/cooler  than  2M0355+11  and  2M2148+40,  but  has  a  nearly  identical  spectrum  to  2M2244+20.  The  Csi  are  well-matched  in 
2M2244+20  and  W0047+68,  but  the  signal-to-noise  is  not  sufficient  to  compare  the  Rb  i  line.  Again,  the  spectra  are  from  different  instruments  so  we  have  smoothed 
the  data  to  match  resolutions. 


but  because  of  the  position  on  the  sky  (/  =  89.0,  b  =  —10.5), 
this  mainly  affects  the  V  component  of  space  motion.  U  is 
strongly  constrained:  if  we  adopt  =  0  ±  40kms_1,  then 
U  =  —34.0  zb  0.5kms_1,  and  W=  —5.6  ±  7kms-1.  This  rules 
out  membership  in  any  of  the  known  very  young  nearby  moving 
groups.  2M2 148+49  is  intrinsically  fainter  at  the  J  band  than 
ordinary  field  L6  dwarfs  and  brighter  at  WISE  bands. 

2.4.  Photometry 

W0047+68  was  observed  as  part  of  normal  survey  operations 
of  the  Pan-STARRS  wide-held  optical/NIR  PS1  survey  (Kaiser 
et  al.  2010).  We  report  AB  magnitudes  for  the  y  and  z  bands, 
which  are  described  in  Tonry  et  al.  (2012).  The  source  was  not 
detected  in  the  shorter  wavelength  bands.  WISE  and  Two  Micron 
Sky  Survey  (2MASS)  photometry  were  previously  reported 
in  Paper  I. 

We  monitored  W0047+68  in  the  J  band  using  the 
FLAMINGOS  camera  (Elston  et  al.  2003)  on  the  Kitt  Peak 
2.1m  telescope  in  2012  October  (NOAO-Program  12B-0233) 
on  three  consecutive  nights  for  four  hours  per  night,  after  the 
main  target  for  the  program  had  set.  We  find  no  evidence  of 
variability  at  the  level  of  0.05  mag  or  greater. 

Although  unobservable  with  WISE ,  many  of  the  directly 
imaged  planets  have  been  observed  at  L'  and  Mf .  W0047+68  has 
not  been  observed  through  these  filters,  but  we  may  estimate  the 
magnitudes  using  WISE  photometry  and  other  L  dwarfs.  Leggett 
et  al.  (2007)  measured  M'  =  11.90  ±  0.03  for  2M2244+20. 


Scaling  by  the  W2  magnitudes,  this  suggests  M '  =  11.04  for 
W0047+68  and  Mw  =  10.6. 

2.5.  Luminosity 

We  can  obtain  the  bolometric  magnitude  of  W0047+68  by  in¬ 
tegrating  the  spectral  energy  distribution,  improving  the  estimate 
from  Paper  I  using  our  new  data.  The  near-infrared  component 
is  well-constrained  by  the  Paper  I  spectrum  normalized  to  agree 
with  the  2MASS  photometry.  In  the  optical,  we  use  the  MMT 
spectrum.  There  is  negligible  contribution  at  less  than  0.6  pm. 
We  use  the  Spitzer  Infrared  Spectrograph  (IRS)  spectrum  of 
2M2244+20  (Stephens  et  al.  2009)  renormalized  by  a  factor  of 
2.23  (from  the  relative  WISE  W1  and  W2  photometry)  in  the 
range  of  5.15  /rm-14.15  pm  and  a  Rayleigh-Jeans  tail  beyond 
that.  This  spectral  energy  distribution  is  shown  in  Figure  6.  The 
wavelength  range  between  2.55  pm  and  5.15  pm  is  constrained 
by  the  WISE  photometry,  but  we  require  higher  resolution  to 
account  for  the  non-smooth  spectrum.  We  consider  the  A-type 
(1200  K)  and  AE-type  (1300-1500  K)  models  of  Madhusudhan 
et  al.  (2011)  and  the  preliminary  BT-Settl  (1300-1500  K)  of 
Allard  et  al.  (2013).  In  each  case,  we  integrate  the  model  spec¬ 
tra  over  the  range  of  2.55-4.00  pm,  normalized  to  the  WISE 
WI  photometry,  and  the  range  of  4.00-5.15  pm,  normalized 
to  the  WISE  W2  photometry.  (Because  the  observed  colors 
are  redder  than  the  models,  we  allow  discontinuities.)  Aver¬ 
aging  the  three  best- fitting  near-infrared  models  (A  1200  K,  AE 
1300  K,  and  BT-Settl  1400  K,  all  logg  =  4.0)  we  obtain  an 
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Wavelength  (microns) 


Figure  6.  Top:  the  observed  spectral  energy  distribution.  The  Pan-STARRS,  2MASS,  and  WISE  photometry  shown  as  solid  circles.  The  solid  lines  show  the  optical 
(MMT)  and  near-infrared  (IRTF)  spectra  of  W0047+68  and  the  Spitzer  IRS  spectrum  of  2M2244  normalized  to  match  W0047+68.  The  dotted  lines  show  a  very  thick 
cloud  (A-type)  model  (Madhusudhan  et  al.  2011)  with  Teff  =  1200 K,  logg  =  4.0,  a  BT-Settl  model  (Allard  et  al.  2013)  with  Teff  =  HOOK,  logg  =  4.0,  and  a 
blackbody  (HOOK)  spectra  normalized  to  match  the  W1  and  W2  photometry;  these  are  allowed  to  have  discontinuities  at  2.55  /z m,  4.00 /zm,  and  5.15  /zm.  Bottom: 
the  cumulative  spectral  energy  fraction.  Half  of  the  total  energy  of  W0047+68  is  emitted  shortward  of  2.9  /zm. 


apparent  bolometric  magnitude  of  16.30  with  an  uncertainty  of 
±0.02  magnitudes  due  to  the  spread  in  models,  all  of  which  fea¬ 
ture  deep  water  bands.  If  blackbodies  without  absorption  bands 
are  used  for  this  region,  the  resulting  bolometric  magnitude  is 
up  to  0.05  mag  brighter,  while  an  observed  AKARI  spectrum  of 
an  ordinary  L8  dwarf  (Sorahana  &  Yamamura  2012;  Sorahana 
et  al.  2013)  normalized  to  the  W1  photometry  is  0.02  mag  fainter; 
other,  higher-gravity  bluer  models  are  within  this  range  of  un¬ 
certainty.  There  is  an  additional  uncertainty  from  the  uncertainty 
in  the  2MASS  and  WISE  photometry  and  other  spectra,  which 
we  estimate  is  ±0.025  mag.  The  derived  W0047+68  luminos¬ 
ity  is  (3.58  ±  0.29)  x  10-5  L0.  The  cumulative  contribution  to 
the  luminosity  is  also  shown  in  Figure  6:  half  of  the  energy  is 
emitted  redward  of  2.9  p m.  Following  the  same  procedure  for 
2M2148+40,  but  with  the  BT-Settl  1500  K  model  yields  an  ap¬ 
parent  bolometric  magnitude  of  15.12  ±  0.03  and  luminosity  of 
L  =  (4.64  ±0.14)  x  10-5Lo. 

3.  PROPERTIES  OF  W0047+68 
3.1.  AB  Dor  Moving  Group  Membership 

The  spectra  of  W0047+68  include  a  number  of  age  indicators. 
The  most  conservative  age  estimate  for  W0047+68  comes  from 
the  fact  that  we  detect  lithium:  for  the  observed  luminosity,  this 
indicates  that  its  age  is  less  than  one  billion  years  old  (Chabrier 
et  al.  2000).  The  near-infrared  spectroscopy,  however,  points  to 
an  even  younger  age,  less  than  that  of  most  lithium  L-type  brown 
dwarfs.  Thompson  et  al.  (2013)  suggested  that  W0047+68  is  a 
“low-gravity  object  and  hence  young,”  based  on  their  analysis 


of  a  moderate-resolution  spectrum  similar  to  the  one  we  have 
obtained.  Our  near-infrared  surface  gravity  classification  of 
INT-G  (Section  2.1)  indicates  an  age  of  ~  100-150  Myr  for 
W0047+68,  but  not  <50  Myr,  though  this  system  is  primarily 
based  upon  warmer  objects.  Many  of  the  youngest  stars  are 
found  to  be  members  of  moving  groups. 

The  three-dimensional  space  velocity  (Figure  7)  is  consistent 
with  the  young  AB  Dor  Moving  Group  (Zuckerman  et  al. 
2004;  Torres  et  al.  2008). 13  Most  ABDMG  members  are  in 
the  southern  hemisphere,  but  in  fact  the  group  is  quite  loose 
and  includes  northern  hemisphere  members.  Given  the  parallax, 
proper  motion,  radial  velocity,  and  lithium  age  constraint,  we 
find  that  the  BANYAN  II  model  (Gagne  et  al.  2013)  assigns  an 
ABDMG  membership  probability  of  99.96%  to  W0047+68.  As 
seen  in  Figure  7,  W0047+68  lies  fairly  close  to  the  group  in 
physical  space,  and  the  peculiar  velocity  necessary  to  move 
W0047+68  from  the  nucleus  of  the  ABDMG  to  its  present 
position  is  <lkms-1.  Furthermore,  as  discussed  in  detail 
in  Section  3.5,  mid-L  dwarf  ABDMG  candidate  members, 
including  the  companion  CD— 35  2722B  (Wahhaj  et  al.  2011) 
which  is  a  definite  member,  are  unusually  red  in  J  —  K  color. 

Therefore,  on  the  basis  of  three  indicators  (three-dimensional 
space  motion,  spectroscopic  surface  gravity,  and  near-infrared 
color),  we  propose  that  W0047+68  is  a  bona  fide  member  of 
the  ABDMG.  One  could  ask  if  it  is  plausible  that  an  ABDMG 
brown  dwarf  member  is  so  close  to  the  Sun.  Torres  et  al.  (2008) 

13  In  Paper  I,  we  noted  that  the  direction  of  the  proper  motion  of  W0047+68 
agreed  with  the  very  young  (~  12  Myr)  ft  Pic  Moving  Group.  The  parallax 
rules  ft  Pic  Moving  Group  membership  out. 
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Figure  7.  Comparison  of  the  position  and  velocities  of  W0047+68  (red  filled  circle)  to  nearby  young  associations  (Torres  et  al.  2008).  Although  at  the  edge  of 
the  AB  Dor  Association  the  kinematic  agreement  is  excellent.  The  size  of  the  W0047+68  symbol  is  the  1  a  uncertainty.  Luhman  et  al.  (2005)  place  the  Pleiades  at 
t/=  —6.6  ±  0.4kms_1,V  =  —27.6  ±  0.3kms_l,  — 14.5  ±  0.3  km  s"1,  close  to  both  the  center  of  the  ABDMG  distribution  and  W0047+68  itself  in  velocity  space, 
but  in  physical  space  the  cluster  is  located  at  X  =  —1 19  pc,  Y  =  28  pc,  and  Z  =  —53  pc. 


list  eight  ABDMG  M  dwarf  members  within  16  pc  of  the  Sun 
and  Kirkpatrick  et  al.  (2012)  find  that  the  ratio  of  brown  dwarfs 
to  stars  is  1 :6,  so  the  presence  of  one  or  two  nearby  ABDMG 
brown  dwarfs  is  not  surprising. 

3.2.  Age  and  Metallicityfor  W0047+68 

We  can  leverage  the  current  knowledge  of  the  ABDMG, 
which  is  based  on  a  wide  variety  of  stellar  masses  and  spectral 
types,  to  understand  the  properties  of  W0047+68  and  the  origin 
of  its  peculiar  spectrum.  Although  early  discussions  of  ABDMG 
(Zuckerman  et  al.  2004)  assumed  an  age  of  ~50  million  years, 
current  evidence  points  to  an  age  of  ~  100-1 25  Myr  (Luhman 
et  al.  2005;  Barenfeld  et  al.  2013).  Indeed,  Luhman  et  al.  (2005) 
argue  that  ABDMG  is  “roughly  co-eval”  with  the  Pleiades 
and  is  physically  related  given  the  similar  space  velocities, 
an  argument  supported  by  Ortega  et  al.  (2007).  This  age  is  in 
excellent  agreement  with  our  spectroscopic  estimate  of  surface 
gravity  intermediate  between  the  youngest  association  and  older 
field  brown  dwarfs.  Of  course,  even  the  age  of  the  rich,  well- 
studied  Pleiades  is  still  somewhat  uncertain  (see  Soderblom 
et  al.  2009),  but  W0047+68’s  age  constraint  is  much  better  than 
that  of  most  known  brown  dwarfs. 

The  relative  importance  of  surface  gravity  and  composition 
is  debated  in  both  unusually  red  brown  dwarfs  and  massive 
gas  giant  planets.  Crucially,  the  composition  of  the  ABDMG 
is  also  known  to  be  near  solar.  Ortega  et  al.  (2007)  compiled 
the  recent  literature  measurement  of  15  ABDMG  members  and 
found  a  mean  [Fe/H]  =  —0.02  ±  0.02.  Viana  Almeida  et  al. 


(2009)  measured  a  mean  [Fe/H]  =  —0.01]  for  12  ABDMG 
members,  with  a  standard  deviation  of  0.09,  but  after  applying 
a  correction  for  a  trend  in  effective  temperature  found  a  mean 
[Fe/H]  =  0.04  with  a  standard  deviation  of  0.05.  They  also 
reported  that  for  their  entire  young  star  sample  the  mean  values 
are  [Si/H]  =  -0.11  and  [Ni/H]  =  0.01.  Biazzo  et  al.  (2012) 
reported  a  mean  [Fe/H]  =  0. 10  db  0.03  for  the  five  stars  in 
the  ABDMG,  and  found  that  10  elements  (Na,  Mg,  Al,  Si,  Ca, 
Ti,  Cr,  Fe,  Ni,  and  Zn)  have  abundance  patterns  that  are  solar 
and  “typical”  of  the  local  Galactic  thin  disk.  Barenfeld  et  al. 
(2013)  claim  a  somewhat  lower  [Fe/H]  =  0.02  ±  0.02  for 
the  ABDMG  nucleus  and  measure  seven  other  elements  (Na, 
Mg,  Al,  Si,  Cr,  Mn,  and  Ni)  to  be  close  to  solar  values.  Since 
these  [Fe/H]  values  are  characteristic  of  the  bulk  of  the  local 
thin  disk,  the  vast  majority  of  field  L  dwarfs  must  have  similar 
values.  Whether  the  scatter  in  values  is  a  matter  of  observational 
uncertainties  or  true  scatter  could  be  debated.  Barenfeld  et  al. 
(2013)  argue  that  they  see  real  scatter  in  composition  in  the 
ABDMG  “stream”  members  and  that,  therefore,  many  of  the 
more  distant  ABDMG  objects  are  not  from  a  common  origin. 
Nevertheless,  even  these  objects  are  all  near  solar  in  [Fe/H]  and 
other  elements,  as  is  the  Pleiades  (Soderblom  et  al.  2009),  so 
we  conclude  that  W0047+68’s  composition  is  constrained  to  be 
near-solar  even  if  the  ABDMG  stream  is  not  truly  homogeneous. 
This  constraint  includes  the  key  elements  involved  in  cloud 
formation:  Fe  and  Al  (along  with  O)  are  the  major  elements  for 
high-temperature  condensates  (i.e.,  in  the  deeper  atmospheric 
layers  of  W0047+68)  and  Mg  and  Si  (and  O)  at  moderate 
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Figure  8.  Color-magnitude  diagram  comparing  W0047+68  and  2M2 148+40 
to  ordinary  field  L  and  T  dwarfs,  very  low  mass  brown  dwarfs,  and  exoplanets 
compiled  by  Liu  et  al.  (2013b). 


temperature  is  7 eff  ~  1270-1300 K  with  log#  ~  4.5.  Burrows 
et  al.  (2011)  report  that  high  metallicity  or  thick  cloud  models 
each  lead  to  an  increase  in  the  radius  over  “standard”  models, 
which  would  lower  the  temperature  below  1400K  even  if 
W0047+68  is  old. 

Given  an  effective  temperature  of  ~1300K,  W0047+68 
shows  that  intermediate  surface  gravity  brown  dwarfs  have  a 
much  different  L/T-type  transition.  For  older,  high  gravity  field 
brown  dwarfs,  this  effective  temperature  is  normally  associated 
with  the  early  T  dwarfs  (Cushing  et  al.  2008;  Stephens  et  al. 
2009)  and  certainly  with  spectral  types  that  have  already  begun 
the  move  to  the  blue  in  J  —  K.  The  early  T  dwarfs  have  methane 
bands  in  the  H  band  and  K  band  (Burgasser  et  al.  2006).  No 
such  features  are  seen  in  W0047+68,  and  we  also  did  not  detect 
the  high  amplitude  7-band  variations  seen  in  ~1300K  early 
T  dwarfs  and  attributed  to  patchy  clouds  (Artigau  et  al.  2009; 
Radigan  et  al.  2012;  Apai  et  al.  2013).  (However,  W0047+68’s 
low  v  sin  i  suggests  that  we  may  be  viewing  it  near  pole-on.) 
We  emphasize  that  this  result  depends  upon  the  evolutionary 
models:  W0047+68  has  a  luminosity  typical  of  that  of  ordinary 
L7  dwarfs  (Vrba  et  al.  2004;  Golimowski  et  al.  2004),  but  its 
young  age  (and  theoretically  larger  radius)  requires  a  lower 
temperature.  As  shown  by  Liu  et  al.  (2013b),  the  younger  object 
PSO  J3 18-22  is  an  even  more  extreme  case,  with  nearly  the 
same  luminosity  and  spectral  type  as  W0047+68,  but  a  deduced 
temperature  of  1160+4°  K. 


temperatures  characteristic  of  the  upper  atmospheric  layers 
(Marley  et  al.  2002;  Helling  et  al.  2008b;  Visscher  et  al.  2010). 

Using  the  Chabrier  et  al.  (2000)  models  and  the  observed 
luminosity,  W0047+68  would  have  logg  &  4.5  and  M  % 
0.018  A/0  for  age  120Myr.  Similar  values  are  found  with  the 
Burrows  et  al.  (1997)  models.  The  low  surface  gravity  and  solar 
metallicity  suggest  that  the  thick  clouds  of  W0047+68  are  a 
normal  feature  of  a  young  brown  dwarf.  In  Figure  8,  we  show 
the  color-magnitude  diagram  for  brown  dwarfs  and  planetary 
mass  objects  (Liu  et  al.  2013b),  updated  with  W0047+68  and 
2M2 148+40.  W0047+68  lies  close  to  PSO  J3 18-22,  adding  to 
the  evidence  of  a  low-gravity  sequence  that  extends  redward  of 
the  normal  L/T  brown  dwarf  sequence. 

3.3.  Effective  Temperature 

In  Paper  I,  we  found  that  although  models  from  different 
groups  agreed  that  unusually  thick  clouds  were  needed  to  fit 
the  near-infrared  spectrum,  the  predicted  effective  temperatures 
ranged  from  as  high  as  1600K  to  as  low  as  HOOK.  With  our 
measurement  of  the  luminosity,  we  can  make  another  estimate 
of  the  effective  temperature  using  theoretical  predictions  of  the 
radius.  We  find 

W-  ® 

The  highest  temperature  and  smallest  radius  for  W0047+68 
can  be  obtained  by  using  the  oldest  evolutionary  models  in  which 
lithium  has  not  been  depleted.  This  would  allow  7eff  ^  1450  K 
and  an  age  of  ~1  billion  years,  with  logg  ^  5.  ABDMG 
membership  implies  a  larger  radius  and  lower  temperature. 
Adopting  an  age  of  125  Myr  for  ABDMG  membership.  Burrows 
et  al.  (1997)  have  Af0  =  0.020Mo  and  R  =  0.117  Re; 
Chabrier  et  al.  (2000)  have  M  =  0.018  M0  and  R  =  0.124  RQ. 
Correspondingly,  our  best  estimate  of  W0047+68’s  effective 


3.4.  Comparison  to  Thick  Cloud  Models 

The  properties  of  condensate  cloud  layers  depend  on  complex 
processes  that  are  difficult  to  model  completely,  including 
nucleation,  growth,  and  settling  (Ackerman  &  Marley  2001; 
Helling  et  al.  2008b).  Each  theoretical  effort  to  predict  synthetic 
spectra  makes  different  choices  in  parameterizing  the  cloud 
properties.  These  were  compared  by  Helling  et  al.  (2008a)  in  test 
cases  aimed  at  modeling  typical  field  brown  dwarfs.  In  Paper  I, 
we  discussed  how  theoretical  models  from  three  different  groups 
(UCM;  Tsuji  2005; /sed,  Ackerman  &  Marley  2001;  Marley 
et  al.  2002;  Stephens  et  al.  2009;  very  thick  (A-type)  and 
intermediate  thickness  (AE-type)  cloud  models,  Madhusudhan 
et  al.  2011)  were  able  to  fit  the  low-resolution  W0047+68 
near-infrared  spectrum.  In  each  case,  models  can  fit  the  J  —  K 
color  and  qualitatively  reproduce  the  spectrum  by  adjusting 
a  free  parameter  to  increase  the  thickness  of  the  condensate 
clouds.  Another  set  of  models,  DRIFT-PHOENIX  (Witte  et  al. 
2009),  was  used  by  Witte  et  al.  (2011)  to  fit  low-resolution 
near-infrared  spectra  of  2M2 148+40  and  2M2244+20.  They 
derived  warm  temperatures  (7ef r  =  1500  K)  and  very  low 
gravities  (log  g  =  3.0  and  3.5,  respectively),  but  discussed  the 
model  grid  limitations  for  such  red  objects  and  concluded  that 
their  temperatures  are  likely  overestimated  and  the  gravities 
underestimated. 

We  can  now  also  compare  to  a  fifth  set  of  models.  Allard 
et  al.  (2012,  2013)  have  presented  a  new  model  family  called 
BT-Settl.  The  mixing  and  cloud  properties  are  informed  by 
two-dimensional  radiation  hydrodynamic  (hereafter  RHD)  sim¬ 
ulations  (Freytag  et  al.  2010).  The  result  is  that  for  cool  L 
dwarfs,  low-gravity  objects  are  predicted  to  be  redder,  with 
thicker  clouds.  We  show  the  W0047+68  observed  spectrum 
with  the  model  predictions  in  Figure  9.  While  our  analysis 
of  W0047+68  favors  7eff  ^  1300K  and  logg  &  4.5,  the 
BT-Settle  grid  requires  an  even  lower  gravity  but  higher  temper¬ 
ature  (7eff  =  1400-1500  K,  log  g  =  4.0)  to  match  the  spectrum. 


7eff  =  (1413K±28  K)| 
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Figure  9.  Comparison  of  observed  W0047+68  spectrum  to  low-gravity  (logg  =  4.0,  4.5)  BT-Settl  models:  each  is  normalized  to  agree  in  the  1.50-1.75  /xm  region. 
While  we  conclude  that  T0 ff  =  1300-1430  K  and  log  g  =  4.5,  the  models  with  log  g  =  4.0  do  a  better  job  of  reproducing  the  red  near-infrared  spectrum. 


(Note  that  the  log  g  =  3.5  BT-Settl  models  are  even  redder  and 
can  also  fit.) 

It  has  long  been  recognized  that  J—  K  colors  are  relatively 
poor  temperature  indicators  for  brown  dwarfs,  even  among  the 
typical  field  population  because  of  the  effects  of  clouds  (Marley 
et  al.  2002;  Knapp  et  al.  2004;  Stephens  et  al.  2009).  Our 
results  again  underline  this  issue.  The  Unified  Cloudy  Model 
and  BT-Settl  give  qualitatively  good  fits  to  W0047+68,  but  are 
warmer  by  100-200  K  than  our  estimate  of  1300  K,  while,  on 
the  other  hand,  the  very  thick  (A-type)  and  /se d  =  1  models 
give  temperatures  of  100-200  K  too  low  when  matching  the 
spectra.  The  plots  in  Figure  9  have  arbitrary  normalizations.  An 
interesting  check  is  to  derive  the  radius  from  the  near-infrared 
region  alone  instead  of  using  our  procedure  in  Section  3.3:  given 
a  model,  the  observed  parallax,  and  the  observed  flux-calibrated 
spectrum,  the  only  adjustable  parameter  is  the  radius.  The 
BT-Settl  results  when  fitting  over  the  range  of  1.2-2. 5  /xm  with 
the  log  g  =  4.0  models  are  0.13  RQ  (1200  K),  0.12  RQ  (1300  K), 
0.94  Rq  (1400  K),  and  0.82  RQ  (1500  K),  consistent  with  the 
discussion  in  Section  3.3.  The  1500  K  model  has  a  qualitatively 
good  fit  but  requires  a  radius  that  is  considerably  smaller  than  we 
expect  for  a  brown  dwarf  with  lithium.  We  believe  the  balance 
of  evidence  is  against  this  possibility,  though  it  should  be  noted 
that  Sorahana  et  al.  (2013)  have  used  their  own  model  fitting 
to  AKARI  data  to  argue  that  (ordinary)  field  late-L  dwarfs  are 
warmer  and  much  smaller  than  predicted.  Cushing  et  al.  (2008) 
showed  that  fits  to  limited  wavelength  regions  can  almost  always 
work,  but  can  give  much  different  effective  temperatures  than 
fits  to  the  entire  spectral  energy  distribution. 


One  important  conclusion  from  our  work  is  that  model  colors 
should  not  be  used  to  derive  surface  gravities  because  they 
cannot  yet  predict  cloud  properties  well  enough.  In  the  cases 
with  a  full  grid  of  surface  gravities  (Madhusudhan  et  al.  2011; 
Allard  et  al.  2012),  one  must  go  to  log  g  =  4.0  or  even  less  to  get 
thick  enough  clouds  to  match  in  the  infrared.  This  would  in  turn 
require  that  W0047+68  be  below  the  deuterium-burning  limit 
and  very  young,  much  younger  than  the  best  estimates  of  the 
ABDMG  age.  For  example,  in  the  Baraffe  et  al.  (2002)  models 
at  20Myr,  a  0.007  M0  would  have  log  g  =  4.0  and  match  the 
observed  luminosity  with  Te ff  =  1200K,  while  a  0.009  M© 
would  have  lo gg  =  4.1  and  reff  =  1400  K  but  would  be  90% 
more  luminous  than  observed.  Manjavacas  et  al.  (2014)  have 
argued  that  the  high  altitude  dust  in  BT-Settl  models  should 
be  increased  to  better  fit  the  warmer  (M9.5-L3)  young  brown 
dwarfs;  the  same  is  evidently  true  for  W0047+68. 

The  optical  spectrum  of  W0047+68  also  poses  a  challenge 
for  our  understanding.  The  spectrum  is  dominated  by  the  broad 
wings  of  the  Ki  which,  in  turn,  are  the  dominant  feature  that 
can  be  used  for  spectral  typing;  these  are  evidently  saturated 
around  1300  K  even  though  their  weakness  is  an  indicator 
of  low  gravity  in  warmer  objects  (Cruz  et  al.  2009).  This 
was  recognized  by  Kirkpatrick  et  al.  (2008)  who  argued  that 
2M2244+20  is  young,  yet  found  a  normal  L  dwarf  spectrum 
(shown  in  our  Figures  3  and  5).  They  speculated  that  “the 
lower  gravity  (optical)  spectrum  will  have  weaker  alkalis  and 
stronger  hydrides  than  a  high-gravity  L  dwarf  of  comparable 
temperature.”  This,  in  turn,  creates  “a  degeneracy  at  mid-  and 
late-L-types  whereby  low-gravity  objects  are  very  difficult  or 


10 


The  Astrophysical  Journal,  799:203  (15pp),  2015  February  1 


0.65  0.70  0.75  0.80  0.85  0.90 

Wavelength  (Microns) 

Figure  10.  Comparison  of  observed  W0047+68  spectrum  to  Madhusudhan  et  al. 
(2011)  intermediate  thickness  clouds  (AE  type)  and  very  thick  cloud  (A  type) 
models  with  low  surface  gravity  (logg  =  4.0)  models.  These  two  models  are 
normalized  to  the  W0047  /7-band  region  (as  in  Paper  I).  The  strength  of  the  broad 
potassium  line  is  under-predicted,  and  the  pseudo-continuum  is  too  high,  due 
to  the  overly  low  surface  gravity.  We  also  show  a  Burrows  et  al.  (2006)  model 
with  1300  K  log  g  =  5.0  and  normal,  thin  (“E-type”)  clouds,  normalized  in  the 
region  0.7-0.9  /zm,  which  demonstrates  that  this  model  family  does  reproduce 
the  main  features  of  L  dwarfs  in  the  optical. 

impossible  to  distinguish  in  the  optical.”  If  our  interpretation 
of  W0047+68  is  correct,  this  must  be  the  case  because  (except 
for  the  strengths  of  the  Cs  and  Rb  lines)  it  is  very  close  to  the 
L7  standard.  Another  important  factor  in  creating  degeneracy 
may  be  additional  opacity  due  to  dust,  which  we  know  must 
be  different  in  the  near-infrared  and  may  be  different  in  the 
optical.  Pavlenko  et  al.  (2000)  discussed  the  optical  spectra  of 
L  dwarfs  and  showed  that  the  broad  K  i  line  and  the  sharper  Cs  i 
and  Rbi  lines  in  the  TQ ff  =  1200-1500  K  range  are  sensitive 
to  surface  gravity  and  dust  extinction;  they  could  change  the 
best-fit  effective  temperature  and  surface  gravity  by  increasing 
or  decreasing  them  together,  and  additional  opacity  helped  the 
quality  of  fits.  Lodders  (1999),  Lodders  &  Fegley  (2002),  and 
Marley  et  al.  (2002)  discuss  the  sensitivity  of  the  Ki  chemistry 
to  sedimentation  and  cloud  formation.  Thus,  it  may  be  possible 
to  find  models  where  the  “L7”  optical  spectra  are  relatively 
insensitive  to  gravity. 

Synthetic  optical  spectra  are  available  for  the  BT-Settl  and 
Madhusudhan  et  al.  (2011)  model  families,  and  have  not  been 
tuned  to  match  observations.  In  Figure  10,  we  compare  the 
predicted  and  observed  far-red  spectra  using  the  Madhusudhan 
et  al.  2011  very-thick  (A)  and  intermediate  thickness  (AE) 
models  normalized  to  agree  in  the  H  band  (as  in  Figure  6 
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of  Paper  I).  In  both  cases,  the  absorption  by  the  broad  Ki 
wings  is  underestimated  and,  as  a  result,  the  continuum  is  too 
high.  Since  the  optical  and  near-infrared  spectra  of  typical  L 
dwarfs  are  fit  well  by  the  “E  type”  thin  cloud  models  and 
log  g  =  5.0  (Burrows  et  al.  2006),  one  would  expect  them 
to  also  fit  the  optical  spectrum  of  W0047+68.  Such  a  fit  is  also 
shown  in  Figure  10.  The  optical  BT-Settl  models  are  shown  in 
Figure  11.  There  are  also  disagreements  with  the  observations, 
with  the  pseudo-continuum  lying  too  low,  but  the  TQ ff  =  1300  K, 
log  g  =  4.0  model  is  reasonably  similar  in  overall  morphology 
and  slope,  though  it  lacks  the  hydrides;  CrH  and  FeH  are  better 
matched  with  the  logg  =  4.5  models.  The  change  in  shape  of 
the  core  of  the  Ki  resonant  doublet  in  the  logg  =  4.0  models 
from  1500  K  to  1300K  is  similar  to  the  differences  between 
2M0355+11  and  W0047+68  (Figure  5),  and  suggests  that  the 
basic  picture  in  which  the  core  of  the  K  i  feature  is  saturated  for 
late  L  dwarfs  even  for  low  gravities  may  be  correct.  In  both  the 
BT-Settl  and  A/AE  models,  the  logg  =  4.0  models  have  no 
CrH  or  FeH  features:  this  prediction  could  be  tested  in  the  ft  Pic 
Moving  Group  planetary  mass  object  PSO  J3 18-22. 

3.5.  AB  Dor  Moving  Group  Brown  Dwarfs 

There  are  now  a  number  of  L-type14  brown  dwarfs  proposed 
as  ABDMG  members.  The  strongest  members,  CD— 35  2722B 
(Wahhaj  et  al.  2011),  2M0355+1 1  (Faherty  et  al.  2013b;  Liu  et  al. 
2013a),  and  W0047+68,  show  a  surprising  spread  in  spectral 
properties.  1RXS  J2351+3127  B  (Bowler  et  al.  2012)  and 
2MASS  01222439  B  (Bowler  et  al.  2013)  still  need  kinematic 
confirmation.  The  similarity  of  2M2244+20  and  W0047+68 
supports  ABDMG  membership  for  2M2244+20;  using  the 
W0047+68  absolute  magnitudes  for  2M2244+20  implies  a 
distance  consistent  with  ABDMG  tangential  motion  (Gagne 
et  al.  2014).  The  signal-to-noise  on  the  2M2244+20  spectrum  is 
too  low  for  a  reliable  measurement  of  lithium.  In  Figure  12, 
we  plot  the  intermediate-resolution  near-infrared  spectra  of 
1RXS  J2351+3127  B,  CD-35  2722B,  2MASS  01222439  B, 
2M0355+11  (Liu  et  al.  2013a),  and  W0047+68. 

Allers  &  Liu  (2013a)  have  discussed  the  discrepancies  be¬ 
tween  2M0355+1 1  (L3  VL-G)  and  CD-35  2722B,  which  is  L3 
INT-G  and  a  companion  to  one  of  the  nucleus  AB  Dor  mem¬ 
bers.  Most  notably,  2M0355+1 1  is  more  strongly  peaked  in  the 
H  band.  Despite  the  very  similar  J  —  Ks  colors  of  W0047+68 
and  2M0355+11,  the  surface  gravity  sensitive  features  are  very 
different,  and  W0047+68  appears  to  be  cooler.  2M0355+11 
is  also  a  magnitude  brighter  in  the  color-magnitude  diagram 
(Figure  8).  The  Pleiades  has  a  similar  age  and  space  motion, 
and  infrared  spectroscopy  shows  early  L  dwarfs  in  the  Pleiades 
are  unusually  red  (Bihain  et  al.  2010).  Allers  &  Liu  (2013b) 
report  that  these  spectra  can  be  classified  as  VL-G  or  INT-G  but 
show  a  spread  in  properties.  Unfortunately,  no  late-type  L  dwarfs 
in  the  Pleiades  are  yet  known,  but  our  adopted  age  of  ABDMG 
implies  that  W0047+68-like  members  should  be  found  in  deep 
Pleiades  proper  motion  searches.  Together,  the  ABDMG  mov¬ 
ing  group  brown  dwarfs  form  an  empirical  intermediate  surface 
gravity  (logg  =4.5)  sequence. 

One  puzzle  is  that  2M0355+11  has  two  very  obvious  indica¬ 
tors  of  low  surface  gravity — the  very  peculiar  optical  spectrum 
(L5 y),  a  very  sharply  peaked  H  band — while  W0047+68  has 
more  subtle  features.  The  inconsistency  of  the  H- band  shape  as 
a  gravity  indicator  is  discussed  in  detail  by  Allers  &  Liu  (2013a). 


14  Delorme  et  al.  (2012)  have  also  proposed  CFBDSIR  2149-0403,  an 
unusually  red  T  dwarf,  as  a  candidate  planetary  mass  ABDMG  member. 
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Figure  11.  Comparison  of  observed  W0047+68  spectrum  to  low-gravity  (log  g  =  4.0,  4.5)  BT-Settl  models:  each  is  normalized  to  agree  in  the  1.5-1.75  /zm  region. 
See  the  discussion  in  Section  3.4. 
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Figure  12.  Comparison  of  ABDMG  L-type  brown  dwarfs.  Note  the  sharply  peaked  2M0355+11  H- band  spectrum  compared  to  the  other  cool  L  dwarf  members. 
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With  regard  to  this  pair  of  objects,  three  possibilities  merit  dis¬ 
cussion.  One  possibility  is  thatboth2M0355+ll  and  W0047+68 
are  normal  ABDMG  members  whose  differences  are  simply  a 
consequence  of  temperature.  Liu  et  al.  (2013b)  have  found  that 
the  luminosity  of  2M0355+11  is  logL/L0  =  —4.23  ±0.11, 
i.e.,  ~60%  more  luminous  than  W0047+68.  This  suggests  that 
it  is  ~170K  warmer  than  W0047+68,  consistent  with  the  two 
spectral  types’  earlier  classification  in  the  optical.  The  optical 
observed  differences  in  the  Ki  core  (Figure  5)  appear  quali¬ 
tatively  similarly  to  the  BT-Settl  offsets  between  the  1500  K/ 
1400  K  models  and  the  1300  K  models  for  logg  =  4.0;  the 
BT-Settl  low-gravity  models  are  most  “peaky”  at  1500  K,  but 
the  //band  becomes  more  rounded  for  lower  temperatures.  The 
A/AE  models  are  also  not  “peaky”  at  the  H  band  for  1300K. 
A  less  likely  possibility  is  both  2M0355+1 1  and  W0047+68  are 
ABDMG  members,  but  2M0355+11  is  actually  a  lower  mass 
object  very  close  to  the  deuterium-burning  limit,  resulting  in 
a  higher  luminosity  and  temperature  than  W0047+68,  but  a 
lower  surface  gravity.  This  kind  of  complex  behavior  near  the 
deuterium-burning  limit  at  ages  up  to  ~  120  Myr  appear  in  all 
published  model  tracks  and  have  recently  been  discussed  again 
theoretically  by  Spiegel  et  al.  (2011)  and  observationally  by 
Allers  et  al.  (2010)  and  Bowler  et  al.  (2013).  Whether  or  not  this 
suggestion  matches  in  detail  depends  on  the  age  of  ABDMG 
and  the  exact  behavior  of  model  tracks  near  the  deuterium¬ 
burning  limit.  The  third  possibility  is  that  either,  or  both,  are 
not  ABDMG  members.  Of  course,  if  2M0355+11  has  a  lower 
gravity  then  it  is  a  younger  interloper,  perhaps  ~ 30-50  Myr  old. 
This  may  not  be  as  unlikely  as  it  seems.  Pr z-Hipparcos  work 
identified  a  “Local  Association”  (also  known  as  the  “Pleiades 
Super-Cluster”  or  “Pleiades  Moving  Group”)  including  many 
young  stars  in  the  age  range  of  ~50-150Myr  (Eggen  1992; 
Jeffries  1995),  with  a  dispersion  of  ±10kms_1.  The  ABDMG 
is  part  of  this  broader  association.  Montes  et  al.  (2001)  iden¬ 
tified  many  young  stars  with  Local  Association  motions,  and 
found  kinematic  “substructure”  as  well  as  a  fairly  large  spread 
in  velocities.  Makarov  (2003)  confirmed  the  Local  Association 
streaming  motion  in  an  X-ray  selected  sample.  Antoja  et  al. 
(2008;  see  their  Figure  13)  find  many  young  (<100 Myr)  stars 
concentrated  in  this  broad  region  of  U-V  space.  More  contro¬ 
versially,  Lopez- Santiago  et  al.  (2006)  argued  that  the  AB  Dor 
group  could  be  divided  into  a  ~ 3 0-50  Myr-old  subgroup  and  a 
^80-120  subgroup.  (There  are  also,  of  course,  old  stars  with 
similar  motion  so  kinematics  alone  cannot  prove  youth.  In  a 
non-age-selected  Hipparcos  sample,  Bovy  et  al.  2009  find  that 
~10%  of  stars  can  be  associated  with  a  Pleiades  moving  group.) 

4.  2M2148+40:  AN  OLDER  DUSTY  L  DWARF 

In  Paper  I,  we  remarked  that  in  our  low-resolution  near- 
infrared  spectrum  W0047+68  resembled  a  slightly  cooler  (or 
one  subclass  later)  version  of  2M2148+40.  Looper  et  al.  (2008) 
discovered  this  brown  dwarf  and  argued  that  spectroscopic  and 
kinematic  evidence  suggests  it  is  not  young,  but  is  rather  an 
older  brown  dwarf  whose  spectral  peculiarities  might  be  ex¬ 
plained  by  high  metallicity.  (However,  they  did  find  2M2 148+40 
has  lithium  with  an  equivalent  width  of  12.1  ±  0.6  A,  see  Fig¬ 
ure  3.)  At  moderate  resolution  significant  differences  from 
W0047+68  are  revealed.  Allers  &  Liu  (2013a)  found  that  its 
moderate-resolution  infrared  spectrum  lacked  signs  of  low  sur¬ 
face  gravity  or  youth,  and  classified  it  as  L6  FLD-G.  This 
spectrum  is  compared  to  W0047+68  in  Figure  1.  With  a  re¬ 
liable,  precise  distance,  we  can  now  revisit  the  kinematic  ar¬ 
guments.  Our  measurement  of  the  17-component  of  its  space 
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motion  (— 34kms_1)  is  not  consistent  with  any  known  young 
moving  groups,  but  is  not  unusual  among  a  sample  of  stars 
that  are  ~0.5-2  billion  years  old  (such  as  dMe  dwarfs  or 
My  <  4  main- sequence  stars;  Reid  et  al.  2002,  their  Figure  15). 
Allende  Prieto  et  al.  (2004)  measured  [Fe/H]  of  all  (118)  stars 
more  luminous  than  My  =  6.5  mag  within  14.5  pc  from  the 
Sun  and  found  that  2%-5%  have  [Fe/H]  >  +0.25.  Their 
most  metal-rich  star  is  HD  125072  which  has  ((/,  V,  W)  = 
(—37.9,  —17.1,  —34.8)  kms-1,  proving  2M2148+40’s  U  mo¬ 
tion  cannot  be  used  to  rule  out  high  metallicity.  Indeed,  the 
U  velocity  is  consistent  with  the  most  common  velocities  of 
the  very  metal  rich  planet-host  stars  in  the  solar  neighborhood, 
which  are  thought  to  be  scattered  from  more  inner  regions  of  the 
Galactic  disk  (Ecuvillon  et  al.  2007).  Santos  et  al.  (2008)  find 
that  the  nearby  star-forming  regions  have  near-solar  metallicity 
in  both  [Fe/H]  and  [Si/Fe],  and  argue  that  high  metallicity  stars 
originate  in  the  inner  galaxy,  and  Wielen  et  al.  (1996)  show 
that  orbital  diffusion  of  metal-rich  stars  explains  the  presence 
of  high  metallicity  stars  locally.  Wielen  (1977)  find  that  the 
typical  change  in  the  Galactic  orbit  radius  will  be  0. 7-0.9  kpc 
in  0.5-1.0Gyr,  the  oldest  ages  consistent  with  the  presence  of 
lithium,  but  presumably  2M2 148+40  would  be  an  example  of 
an  unusually  large  change. 

Given  the  observed  luminosity,  Te ff  =  1507  ±  11  K  for 
R  =  0.1  Rq.  Applying  the  DUSTY  (Chabrier  et  al.  2000)  evo¬ 
lutionary  models  to  this  luminosity  and  the  observed  presence 
of  lithium,  2M2148+40  would  have  M  ~  0.04-0.055  M0  and 
log g  =  5. 0-5. 2  for  ages  0.5  to  1.0 Gyr.  Our  luminosity  for 
2M2 148+40  is  similar  to  those  of  “normal”  field  L6  and  L7 
dwarfs  (Vrba  et  al.  2004;  Golimowski  et  al.  2004),  so  we  nec¬ 
essarily  derive  similar  effective  temperatures.  In  summary,  the 
new  astrometry  is  consistent  with  the  existing  paradigm  that 
2M2 148+40  has  a  normal  surface  gravity  and  effective  tem¬ 
perature  for  a  lithium  L6-type  dwarf,  so  that  high  metallicity 
(or  else  some  other  unknown  factor)  is  needed  to  explain  its 
unusually  thick  condensate  clouds. 

5.  CONCLUSIONS 

Our  new  measurements  of  W0047+68  confirm  the  view  that  it 
is  an  intermediate  gravity  (log  g  ~  4.5),  young  brown  dwarf  in 
the  ABDMG.  Relatively  bright  with  a  reliable  age,  metallicity, 
and  luminosity,  it  is  a  key  benchmark  for  characterizing  the 
effects  of  clouds  in  substellar  atmospheres.  No  theoretical  model 
is  fully  satisfactory  in  matching  both  the  optical  and  near- 
infrared  spectra  and  luminosity,  but  many  features  are  explained 
and  the  need  for  increased  mineral  clouds  is  clear.  The  low 
effective  temperature  (M300  K)  of  W0047+68  reinforces  the 
finding  that  the  standard  L/T  transition  sequence  for  field  brown 
dwarfs  (logg  ~  5)  does  not  apply  at  intermediate  gravities 
(log  g  =  4.5),  just  as  is  already  known  for  the  low  gravities 
of  planetary  mass  objects  (logg  ~  4.0),  and  the  near-solar 
metallicity  of  the  ABDMG  demonstrates  that  gravity  alone  is 
responsible  for  the  differing  cloud  properties.  The  luminosity 
and  temperature  of  the  older  L  dwarf  2M2 148+40  is  warmer  than 
the  L/T  transition,  but  another  factor,  presumably  metallicity, 
places  it  in  the  redder  sequence  of  the  color-magnitude  diagram 
among  the  young  objects. 

We  are  grateful  to  the  USNO  CCD  parallax  team  for  allowing 
us  to  use  the  preliminary  parallax  measurements.  We  thank  Eric 
Mamajek  for  discussions  of  local  associations;  Kimberly  Aller 
for  examining  Pan-STARRS  data;  Denise  Stephens,  Satoko 
Sorohana,  and  Dagny  Looper  for  providing  electronic  copies 
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